Fluctuation-dissipation relation in a sheared fluid 
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In a fluid out of equilibrium, the fluctuation dissipation theorem (FDT) is usually violated. Using 
molecular dynamics simulations, we study in detail the relationship between correlation and response 
functions in a fluid driven into a stationary non-equilibrium state. Both the high temperature fluid 
state and the low temperature glassy state are investigated. In the glassy state, the violation of 
the FDT is similar to the one observed previously in an aging system in the absence of external 
drive. In the fluid state, violations of the FDT appear only when the fluid is driven beyond the 
linear response regime, and are then similar to those observed in the glassy state. These results are 
consistent with the picture obtained earlier from theoretical studies of driven mean-fleld disordered 
models, confirming the similarity between these models and simple glassy systems. 
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In the past years, a large body of theoretical 0], ex- 
perimental and numerical work has been de- 
voted to the study of correlation and response functions 
in non-equilibrium glassy systems. Generally, the non- 
equilibrium situation of interest is generated by a quench 
below the glass transition temperature. The system falls 
then out of equilibrium, in the sense that its relaxation 
time becomes far greater than the experimental time 
scale. Its properties therefore depend on the time 
elapsed after the quench, also called waiting time. The 
most spectacular dependency is observed for the two-time 
correlation functions which depend then upon the two 
time arguments and not on the time difference only, as 
in an equilibrium system. These correlations decay with 
a relaxation time tr that increases with (aging), typ- 
ically tr (X The same aging behavior is observed in 
the response functions, and a useful quantification of the 
non-equilibrium behavior is encoded in the way the usual 
equilibrium fluctuation dissipation theorem (FDT) is vi- 
olated j|] . For atomic systems of the type studied in this 
paper, this violation has been shown by numerical stud- 
ies to be similar to the one observed in high dimensional 
disordered mean- field models [Q. 

These studies, however, focus on the case of a non- 
equilibrium situation following a quench and on the sub- 
sequent aging phenomena. In Ref. Cugliandolo et 
al suggested another approach to non-equilibrium sys- 
tems, in which the non-equilibrium state was generated 
by "stirring" the system. In such a situation, in which 
energy is fed into the system at a constant rate, the the- 
ory predicts that a stationary non-equilibrium state is 
reached, even when the unstirred system is in a glassy 
state: aging is stopped Further theoretical stud- 

ies on driven mean-field disordered systems estab- 
lished a detailed picture of the non-equilibrium behav- 
ior, which can be summarized as follows. The relax- 
ation of these systems is a two-step process composed 
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of a fast part which is essentially unaffected by the driv- 
ing force, and a slow relaxation occurring on a time scale 
which is a decreasing function of the drive intensity. Si- 
multaneously, a two-temperature pattern appears, the 
fast modes being equilibrated at the bath temperature 
T, while the slow ones have an effective temperature 
Tcff > T. Quantitatively, this is shown by studying a cor- 
relation function C(t) — {A{to+t)B{to)) between two ob- 
servables A and B, and the associated response function 
x{t) = S{A{to+t))/Sh{to), where h is the field conjugated 
to B. Whereas at equilibrium, these quantities are re- 
lated by the FDT, dtC{t) = -ksTxit), this FDT has to 
be generalized in the driven system by introducing an ef- 
fective temperature T^s, through dtC{t) = — fcBTeffx(t). 
In the limit of zero drive, one finds Tcff — > T if T > Tc, 
{Tc is the temperature at which the relaxation time of the 
undriven system diverges) and equilibrium properties are 
recovered. For T < Tc, the limiting effective temperature 
coincides with that of the system aging at the same tem- 
perature [§,0. In Ref. dl it was shown that, for mean 
field models, Tcff is related to the configurational entropy 
available to the system near one free energy minimum, 
confirming the interpretation of Toff as a true temper- 
ature [^,^, both in the thermodynamic and dynamical 
sense. 

In this paper, we investigate the non-equilibrium situa- 
tion created by a steady, homogeneous shear imposed on 
a simple glass forming liquid, i.e. a simple realization of 
the stirring systems considered in Ref. Q which was sug- 
gested by Ref. 0. Note however that, in practice, such a 
situation is more easily realized in a "soft" glassy system 
(a complex fiuid with glassy behavior |l^), which can 
support homogeneous shear flow, than in a usual molec- 
ular glass, in which shear banding and fracture tend to 
take place. Simulations, on the other hand, allow to cre- 
ate an homogeneous shear even in a simple system. The 
goal of this paper is then to determine the behavior of a 
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simple glassy system subjected to a shear and to compare 
the results with the main predictions of mean-field cal- 
culations, keeping in mind the possible relevance to soft 
glassy materials. Although the mean-field scenario is not 
a priori expected to apply to a three dimensional fluid, it 
has the advantage of providing precise predictions, which 
can be compared easily to experimental or numerical re- 
sults. Our focus will be on the fluctuation-dissipation 
relation, the general rheological properties already dis- 
cussed by Yamamoto and Onuki |l^] will be only briefly 
considered. 

The system simulated in this work is a 80:20 mixture of 
2916 Lennard-Jones particles, with interaction parame- 
ters that prevent crystallization |^2|. In all the paper, the 
length, energy and time units are the standard Lennard- 
Jones units (TAA (particle diameter), eaa (interaction en- 
ergy), and To = (mAcri^/e^yi)^/^ [Hi, where tua is the 
particle mass and the subscript A refers to the major- 
ity species . The system has been described in detail 
elsewhere, and its equilibrium (high temperature) prop- 
erties have been fully characterized. At the reduced den- 
sity p = 1.2, where all our simulations are carried out, 
a "computer glass transition" is found in the vicinity of 
Tc ~ 0.435 and the slowing down of the dynamics seems 
to be described well by Mode-Coupling theory The 
aging behavior of the system below this temperature has 
also been characterized extensively, including the viola- 
tion of the FDT in the glassy phase 

The homogeneous shear state corresponding to a pla- 
nar Couette flow is obtained by using the SLLOD algo- 
rithm supplemented by Lees-Edwards boundary condi- 
tions jl5|. The velocity gradient is in the z direction, 
and the fluid velocity in the x direction. Constant tem- 
perature conditions are ensured by thermostatting the 
velocities in the direction perpendicular to the flow using 
the Nose- Hoover method 1 15[ , or in some cases through a 
simple velocity rescaling. The shear rate, denoted by 7, 
naturally introduces a new time scale 7"^ into the prob- 
lem. Obviously, a simulation involving a steady shear 
state is possible only if the available simulation time is 
significantly larger than 7"^. This limits our study to 
shear rates larger than typically 10~"*to, corresponding 
to 10^ time steps. 

The first important consequence of applying a shear 
to the system is that, in accordance with theoretical ex- 
pectations, a non-equilibrium stationary state is reached 
after a transient of a few 7~^. This is true even at tem- 
peratures at which the corresponding undriven system 
behaves like a glass and does not reach equilibrium. The 
benefit is that many of the difficulties associated with the 
simulation of glassy systems, such as the dependence on 
the preparation method of the sample (especially cool- 
ing rate), are eliminated. At a given temperature, a 
sheared sample can be prepared either by cooling an- 
other sheared sample from a higher temperature, or by 
shearing a sample quenched at zero shear. The results 



will be identical provided the first steps of the simula- 
tion are discarded. That the time translation invariance 
property is recovered under shear is illustrated in Fig. 
1, which displays the incoherent scattering function of a 
quenched and sheared sample for several different waiting 
times after the quench. As usual, the incoherent scatter- 
ing function for the particles of type A and wave-vector 
q is defined by 
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In this study, only values of q in the y direction (i.e. per- 
pendicular to both the velocity and the velocity gradient) 
will be considered. This choice avoids the complications 
due to convection by the average flow that arise for other 
wave- vectors IrLll. 




FIG. 1. Incoherent scattering functions for the A particles 
at T = 0.3 and q = 7.47ey, the location of the main peak in 
the structure factor. Bold curves, from left to right: 7 = 0.01, 
0.001 and 0.0001. For 7 = 0.001 five different two-time corre- 
lation functions, taken with time origins equally spaced dur- 
ing a run of duration IO^'tq are also shown as light continuous 
curves. The absence of aging is illustrated by the fact that 
these curves coincide with the bold one. 



In Fig. 0, it is clearly seen that the relaxation time 
tr of the correlations is shear rate dependent. As al- 
ready shown in Ref. [pl| , the viscosity rj scales roughly 
as tr- Consequently, decreases when 7 increases: this 
is a shear-thinning behavior. The viscosity is defined by 
77(7) = (Jxz{i)/i^ where axz is the off-diagonal compo- 
nent of the stress tensor, and is shown in Fig. ||. The 
same type of rheological behavior was obtained in Ref. 
| pT[ for a similar system. A Newtonian regime, where r] 
is independent of 7, is obtained when 7^^ ^ t^. For 
T < 0.45, no such regime is observable. The shear- 
thinning behavior is well characterized by a power law 
^(7) ^ 7~"'"^'. For the lowest temperature investigated 
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here, T = 0.3, one finds a ~ 0.9. This shear-thinning ex- 
ponent is slightly smaller than the one obtained in 
There is however no reason to expect a universal value 
for a, which experimentally is found to be system depen- 
dent, with reported values between 0.5 and 1 [ p^ . 




0.0001 0.001 ■ 0.01 

7 

FIG. 2. Viscosity versus shear rate for temperatures (from 
bottom to top) T = 0.5, 0.45, 0.4 and 0.3. The solid lines are 
guides to the eye, the dashed line corresponds to a power law 



The fact that the time translation invariance is reestab- 
lished by shearing the system implies a superficial resem- 
blance with equilibrium systems. The correlation func- 
tion displayed in Fig. 0, for example, is very similar 
to the correlation functions in a fluid system slightly 
above the glass transition. In the following, the dif- 
ferences between the dynamics of fluctuations in the 
driven and equilibrium system is studied through the 
fluctuation-dissipation relation introduced in the intro- 
duction. In the present case, the two observables un- 
der study are Aq{t) — l/TV^ exp[iq • rj{t)], and 
Bq{t) = J2j cos[q • TCj{t)], where ej is a random vari- 
able taken from a bimodal distribution e^- = ±1. It is 
straightforward to show that after averaging over the dis- 
tribution {cj}, the correlation function (^q(i+to)i?q(to)) 
is equal to Cq(t)/2, where Cq is the incoherent scatter- 
ing function defined in Eq. (^. To compute the re- 
sponse function, a term AH — —hBqit) is added to 
the Hamiltonian. The response function is then Xq(^) — 
5{Ac^{t + tQ))/Sh{to). The procedure to study the FDT 
is then the following. The system is made stationary at 
a fixed temperature and shear rate. The field is switched 
on at to and the observable Aq{t + to) is monitored. The 
same procedure is repeated for several (20 to 80) realiza- 
tions of the charge distribution. This gives the integrated 
response function A/q(i), defined as: 

Mq(t)=^'^'"dt'Xq(t')-^^- (2) 



The last equality holds in the linear response regime, that 
we carefully checked, by choosing a small amplitude for 
the field h (between 0.05 and 0.2). This procedure was 
carried out for four different values of the temperature 
(T = 0.5, 0.45, 0.4 and 0.3) both above and below the 
computer glass transition temperature. At each temper- 
ature, shear rates from 5.10"^ to 10~^ were considered. 
The wave vector was q = 7.47ey. The results are most 
easily analyzed by considering parametric plots of Mq(t) 
versus Cq(t). The slope of these curves is, by definition. 
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FIG. 3. Parametric plots for (a) T = 0.3 and (b) T = 0.5, 
and various shear rates. In both figures, the dashed line is 
the FDT, and has a slope —1/T. The full lines are linear fits 
to the data for 7 > 0. 

In Fig. ^, such plots are shown for (a) T = 0.3, which 
is deep in the glassy region, and (b) T = 0.5, where the 
system is at equilibrium in the absence of shear. For 
T = 0.3, a very clear deviation from FDT is observed 
when the correlation falls below C ~ 0.8 = q, which 
corresponds to the plateau value in the correlation func- 
tions (see Fig. The parametric curve can be very 
well approximated by two straight lines, one with slope 
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— l/ksT (FDT) at large correlations and one with a slope 

— l/fcsTeff for C < q. The latter slope saturates for the 
smallest shear rates to the value T^s — 0.62 > Tc (this 
straight line correctly fits the two values 7 = 0.001 and 
0.0001). Remarkably, this value for Toff is very close (cer- 
tainly within error bars) to the one obtained in Ref . , 
when studying the aging system. Such a coincidence was 
expected on theoretical grounds [^,0 , and exemplifies the 
deep meaning of the effective temperature. As it is in 
fact a property of the free energy landscape available to 
the system, it naturally arises in both types of dynamics 
(aging and driven) of the system. The two-temperature 
pattern is also consistent with earlier expectations that a 
simple atomic glass has the same dynamical behavior of 
mean-field disordered models of the "p-spin" type (p > 2) 
[0,0 . It is worth noting that the value q below which the 
violation of FDT is observed is, again as expected from 
theory, approximately equal to the plateau value of the 
correlation functions. This had not been found in Ref. 

, and the present results show that this was probably 
due to preasymptotic effects in the aging results. 

At T = 0.5, the FDT holds at zero shear (this is 
very well verified numerically, see Fig. ^) and remains 
valid for small enough shear rates. Deviations are ob- 
served only for shear rates large enough to induce non- 
Newtonian behavior and shear-thinning. In this regime, 
it is natural to expect that deviations to FDT will be- 
come more important for increasing shear rates. The 
data clearly indicate that the parametric plot can still be 
fitted by two straight lines, with a shear rate dependent 
slope for the non-FDT part. If this slope is interpreted 
as an effective temperature, this temperature increases 
with increasing the strength of the driving force, which 
is a quite intuitive effect. Again, this is very similar to 
the results obtained in ^ for the disordered mean-field 
model. Finally, we note that the results obtained for 
T = 0.4 and T = 0.45 are similar to those obtained for 
T = 0.3 and T = 0.5, respectively. 

In this paper, we have presented the first numerical 
study of the fluctuation-dissipation relation in a realis- 
tic model of a sheared fluid. This study demonstrates 
that the non-equilibrium fluctuations obey the two-time 
scales, two-temperature scenario which was previously 
derived for mean-field glass models 0] . Although the sys- 
tem is a three dimensional one, which will certainly differ 
from mean-field systems in many respects, the mean-field 
approach provides here a useful conceptual framework to 
rationalize the numerical findings. On short time scales, 
the fluctuation-dissipation temperature is equal to the 
microscopic temperature. As the driven system is forced 
to explore phase space further (slow relaxation) , the cor- 
relation drops below its plateau value and a different 
fluctuation-dissipation temperature is observed, which is 
only weakly dependent on the bath temperature 
This last time scale is strongly shear rate dependent, as 
is well known from rheological experiments. 



Our results were obtained on relatively short time 
scales, for a system which is usually accepted as a rea- 
sonable microscopic equivalent of molecular or metallic 
glasses. We believe, however, that they would be most 
easily applicable and experimentally testable for complex 
fluids, in which the glassy state can be more easily dis- 
rupted to establish a permanent shear flow. 
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